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We demonstrate using computer simulations that the simplest vortex ratchet system for type-II 
superconductors with artificial pinning arrays, a simple asymmetric potential array, exhibits the 
same features as more complicated two-dimensional (2D) vortex ratchets that have been studied in 
recent experiments. We show that the simple geometry, originally proposed by Lee et al. [Nature 
400, 337 (1999)], undergoes multiple reversals in the sign of the ratchet effect as a function of 
vortex density, substrate strength, and ac drive amplitude, and that the sign of the ratchet effect 
is related to the type of vortex lattice structure present. Thus, although the ratchet geometry has 
the appearance of being effectively ID, the behavior of the ratchet is affected by the 2D structure 
of the vortex configuration. When the vortex lattice is highly ordered, an ordinary vortex ratchet 
effect occurs which is similar to the response of an isolated particle in the same ratchet geometry. 
In regimes where the vortices form a smectic or disordered phase, the vortex- vortex interactions are 
relevant and we show with force balance arguments that the ratchet effect can reverse in sign. The 
dc response of this system features a reversible diode effect and a variety of vortex states including 
triangular, smectic, disordered and square. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 

When an overdamped particle is placed in an asymmet- 
ric potential and an additional ac drive is applied, a net 
dc drift velocity or rectification can occur which is known 
as the ratchet effect. Stochastic ratchets can be con- 
structed with Brownian particles, while ordinary ratch- 
ets can be created in deterministic systems 1]. Typically, 
an applied ac drive or periodic flashing of the potential 
couples with some form of asymmetry in the substrate, 
breaking the symmetry of the particle motion. Ratchet 
effects have been studied in the context of molecular mo- 
tors [l| as well as in the motion of colloidal particles Q , 
granular materials ^ and cold atoms in optical trap ar- 
rays 4^. 

Vortices in type-II superconductors act as overdamped 
particles moving in nanostructured potential landscapes. 
The possibility of realizing a ratchet effect to transport 
superconducting vortices was originally proposed by Lee 
et al. who studied vortices interacting with a sim- 
ple asymmetric or sawtooth substrate potential and an 
additional external ac drive. This type of asymmetric 
potential could be fabricated by creating a superconduct- 
ing sample that has an asymmetric thickness modulation. 
The external ac drive is provided by either an oscillating 
external magnetic field or by an applied ac current which 
produces an oscillating Lorentz force on the vortices. The 
numerical simulations of Ref. Q indicated that a net dc 
motion of vortices could arise upon application of an ac 
drive, and that the net vortex drift is in the easy direc- 
tion of the asymmetric substrate, as would be expected 
for a single particle interacting with such a potential. The 



original work on vortex ratchets was followed by a series 
of proposals for alternative vortex ratchets based on 2D 
asymmetric pinning array geometries, each of which pro- 
duces similar types of vortex ratchet effects with the net 
vortex flow in the easy direction of the underlying poten- 
tial [fi,iii[ll. These ratchets operate either through 
a bottleneck effect @, |^ or by means of an effectively ID 
mechanism 0, H, [l3l- Similar ratchet effects have also 
been found for superconductors with asymmetric surface 
barriers (Tlj . 

Recently, new types of 2D periodic pinning geometries 
have been proposed that exhibit not only the regular 
ratchet effect but also a reverse ratchet effect, in which 
the net dc vortex motion is in the hard substrate direction 
[H, [H, [13, [H, some cases, ratchet reversals 

occur when the vortex-vortex interactions at higher fill- 
ings become relevant, resulting in an effective asymmetric 
potential oriented in the direction opposite to the pinning 
structure. As a result, some of the vortices move in the 
hard direction of the pinning substrate, and the net mo- 
tion is described by either an effectively ID [l3)[lM[I3 or 
2D mechanism. Other vortex ratchet systems have 
multiple sign reversals due to symmetry breaking from 
multiple circular ac drives or three state pinning po- 
tentials [l^. The reverse vortex ratchet effect has been 
observed in experiments with 2D arrays of triangular pin- 
ning sites [Ts'l and other types of 2D arrays 'l?]. In the 
work of Silva et al. 17], the periodic pinning array in 
the experimental sample has a two component unit cell 
containing a large pinning site with a small pinning site 
placed to one side. This system exhibits a remarkable 
field dependent multiple reverse ratchet effect. For the 
first matching field and higher odd matching fields, a 
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FIG. 1: Schematic of the simulation geometry. Vortices are 
subjected to an applied current J = Jy which produces a 
driving force F*^^' — Fdc^ or FacX, illustrated in the top por- 
tion of the figure. The bottom portion of the figure shows 
the periodic potential U{x) through which the vortices move. 
The width of the easy direction of the pinning potential is Z"*" 
and the width of the hard direction is Z~, where > l~ and 
+ l~ = a, the pinning lattice constant. 
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FIG. 2: The dc non-normalized average velocity {Vx)Nv 
vs Fdc for varied and both directions of dc drive. The 
curves for Ff^' — xFdc have {Vx)Nv > 0, while the curves 
for Ff"" = -xFdc have {Vx)Nv < 0. In order of increas- 
ing \{Vx)N^\, the curves have n„ = 1.04/A^ 1.22/A^ 1.39/A^ 
1.56/A2, 1.74/A^ 1.91/A^ 2.08/A^ 2.43/A^ 2.78/A^ 3.13/A^ 
3.47/A^ and 3.82/A^. For n„ < 2.08/A^, the negative drive 
critical depinning force is larger than the positive drive crit- 
ical depinning force, > f^, but for n„ > 2.08/A^, this 
reverses and /~ < ■ 



regular ratchet effect occurs, but at even matching fields 
the ratchet effect changes sign. Simulations indicate that 
the pinning sites are capturing multiple vortices, and that 
this can change the direction of the asymmetric effec- 
tive potential at every matching field, as described by a 
ID model. A ratchet effect has also been observed in 
a Josephson-junction array structure which has spatial 



asymmetry [18|. Here, a diode effect occurs where the 



depinning field is higher in the hard direction. A ratchet 
and reversed ratchet effect are also seen for matching 
fields less than one. 

In this work, we perform a detailed investigation of the 
simplest vortex ratchet system created by an asymmetric 
periodic modulation, which was originally shown to ex- 
hibit only an ordinary ratchet effect By considering 
magnetic fields much higher than those used in Ref. Q, 
where the effectively single vortex regime was studied, 
we show that this system can also exhibit ratchet sign 
reversals as a function of vortex density when the vor- 
tex interactions become relevant. Unlike the majority of 
previously studied vortex ratchets, the behavior of this 
system is determined by the detailed 2D structure of the 
vortex configuration, and cannot be explained in terms of 
a simple ID model. We find that the dc response exhibits 
a reversible diode effect correlated with the regions where 
the ratchet effect reverses. The sign of the ratchet effect 
is shown to be related to the nature of the overall vortex 
lattice structure. At magnetic field densities where the 



vortex lattice has a more ordered or crystalline structure, 
the effective vortex interactions cancel or are reduced by 
the symmetry of the lattice, and the response is the or- 
dinary ratchet effect as in the case of a single vortex. 
At fields where the vortex lattice is disordered or has an 
intrinsic asymmetry, the interactions between individual 
vortices become relevant and a portion of the vortices 
move in the easy direction of the pinning potential, re- 
sulting in a reversal of the sign of the ratchet effect which 
can be predicted from force balance arguments. In addi- 
tion to the ratchet response of this system, we show that 
a rich variety of distinct vortex phases can be realized as 
a function of density, including crystalline, smectic, dis- 
ordered and square phases. Commensuration effects also 
create oscillations in the critical depinning force which 
are distinct from those observed for two dimensional pe- 
riodic pinning arrays. 



II. SIMULATION 

We consider a 2D system of size — Ly ^ 24A, where 
A is the London penetration depth, with periodic bound- 
ary conditions in the x and y directions. The sample 
contains Ny vortices interacting with an asymmetric sub- 
strate U (x) of period a = A, illustrated schematically in 
Fig. [T] The vortex density Uy = Ny/L^Ly. A given 
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vortex i obeys the overdamped equation of motion 



V- 



dt 



(1) 



The damping constant rj = 0QC?/27r^^pAr, where d is the 
sample thickness, ^ is the coherence length, is the 
normal-state resistivity, and (po — h/2e is the elementary 
flux quantum. The vortex-vortex interaction force is 



(2) 



j)/rij, and 



Here = |r^ - r^l, f,y = 

the position of vortex i (j). Force is measured in units 
of /o = (j)Q/2Tr^oX^ , and time in units of t = ?7//o- 
Ki{rij/X) is the modified Bessel function which falls off 
exponentially for > A. For computational efficiency, 
the vortex-vortex interaction force is cut off at 6A. We 
have previously found that using longer cutoff lengths 
produces negligible effects [1^. The pinning potential is 
modeled as a sawtooth as shown in Fig. [U 
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Here, Ap — l.O/o and Xi = • x. The width of the long 
or "easy" side of the pinning potential is = (2/3)a 
and the width of the short or "hard" side is l~ = (l/3)a 
such that 1+ + r = a. The term Ff""* in Eq. 1 is the 
force from a dc or ac external driving current. For the 
case of an applied dc drive, Ff^* = F^cX, where Fdc is 
increased from zero in increments of SF^^ — 0-01 every 
10'* simulation time steps. For the case of an applied ac 
drive, F^^* = ±Fac^, where the positive sign is used dur- 
ing the first half of each period r and the negative sign is 
used during the second half period, resulting in a square 
wave centered about zero. The initial vortex configura- 
tions are obtained through simulated annealing, where 
we set F^^* = and apply a thermal force with the prop- 
erties (Ff ) = and {¥f {t)¥j{t')) = 2r)kBT6^j5{t - t'). 
The simulated annealing begins at a high temperature 
F'^ = 5.0 which is well above the melting tempera- 
ture of the vortex lattice, and then F^ is slowly re- 
duced to zero. After annealing, we set F^ = 0, ap- 
ply a dc or ac drive, and measure the average velocity 
(Vx) = {{N~^)Y,f^f)Vi ■ x). For dc driving, we identify 
the critical currents in the positive and negative driving 
directions, and /~, which we define as the drive at 
which 1(14)1 =0.1/o. 



III. DC DEPINNING AND VORTEX LATTICE 
STRUCTURES 

We first compare the dc response of the system for 
driving in the positive (easy) and negative (hard) x di- 
rection. The top set of curves in Fig. [5] show {Vx)Ny 
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FIG. 3: The depinning force curves vs n„ for the system in 
Fig. 2 for driving in the negative direction, (black circles), 
and in the positive direction, (open circles). The diode 
effect reverses near n„ = l.Ql/A^ and n„ = 9.90/A^. Several 
oscillations in the depinning force also occur. 



versus Fdc for F^^* = iiFdc for samples with Ny rang- 
ing from Uy = 1.04/A2 to 3.82/\'^. The bottom set of 
curves show {Vx)Ny versus Fdc for the same samples 
with F^^* — —xFdc- A clear diode effect occurs for 
Uy < 2.08/A^ when the negative (hard) critical depinning 
force is larger than the positive (easy) critical depin- 
ning force /+. Here, fc /fo ^ Ij close to the value of the 
pinning force in the hard direction, while f^/fo ~ 0.5, 
nearly the same as the easy direction pinning force. Thus, 
for Uy < 2.08/A^, the depinning force in each direction 
is primarily determined by the pinning substrate and the 
vortex- vortex interactions are only weakly relevant. A 
pronounced drop in from fc/fo ~ 1 to fc/fa ~ 0.4 
occurs above Uy = 2.08/A^. For Uy > 2.08/A^, the diode 
effect is reversed since now < f^. 

When Uy = 2.08/A^, the negative depinning curve in 
Fig. [2] shows a clear two step depinning process with an 
initial depinning near Fdc — 0.45/o followed by a second 
jump in {Vx) near Fdc — 0-8 fo- At the low first depinning 
threshold, only a portion of the vortices depin, while at 
the second threshold near Fdc = O.8/0, the entire lattice 
depins. 

The values of and /~ undergo changes as a func- 
tion of Uy that are difficult to discern in Fig. [21 In Fig. [3] 
we plot both /+ and versus n„ for systems with up 
to 7000 vortices and Uy — 12.15/A^. There is a posi- 
tive diode effect for «.„ < 1.91/A^ when /~ > /+. The 
diode effect reverses for 1.91/A^ < ny < 9.90/A^, where 
fc < fc'- A second reversal occurs for n„ > 9.90/A^, 
when drops back below f^. The high n„ rever- 
sal is shown in more detail in Fig. 31 where we plot 
(Vx) versus Fdc for two different samples. In Fig. HJa), 
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FIG. 4: (T4) versus Fdc curves for positive driving, F^^* = 
^Fdc (upper dashed line) , and negative driving, F°^* — —5iFdc 
(lower solid line), (a) A sample with n„ — 8.68/A^ showing 
a negative diode effect with /~ < f^. (b) A sample with 
rit, = 10.42/A^ showing a positive diode effect with > ft ■ 



'dt Uy = 8.68/A'^, < /+ and we find a negative 
diode effect. The situation is reversed in Fig. [^b) at 
Uy = 10.42/A^, where /~ > and a positive diode 
effect occurs. 

Fig. [3] indicates that /+ and /~ both undergo oscil- 
lations. A simple ID picture for the vortex structure 
would predict maxima to appear in the depinning force 
at the commensurate fields at which an integer number 
of vortex rows fit inside each pinning trough, given by 
ny = V? j}?, where n is an integer. We do observe a 
maximum in both and f^ near Uy — 4.0/A^, corre- 
sponding to the presence of two ordered commensurate 
vortex rows per channel. Although has a very broad 
plateau centered around Uy — 9.0/A^ which could be in- 
terpreted as the n = 3 matching field, the inadequacy of 
a simple ID picture is indicated by the fact that for /+, 
the corresponding maximum falls at the much lower field 
Uy = 7.81/A^. For vortex densities just above and below 
each of the commensurate fields, the addition or subtrac- 
tion of a small number of vortices does not change the 



number of rows in each trough but instead causes a disor- 
dering of the rows relative to one another, resulting in a 
smectic structure. This occurs because the shear modu- 
lus is smaller than the compression modulus of the rows. 
As larger numbers of vortices are added, the compression 
of each row grows until the rows begin to buckle. The 
buckling transition destroys the alignment perpendicu- 
lar to the vortex rows, causing the vortex configuration 
to become more isotropically disordered. This is similar 
to the commensuration effects found for vortex lattices 
in multilayer systems [2l[. In this case, the field is ap- 
plied along the direction of the layers so that the vortices 
effectively interact with a potential that is periodically 
modulated in one dimension. These experiments found 
oscillations in the magnetization, which is proportional to 
the depinning force. The oscillations were also attributed 
to transitions in the vortex lattice structure from n chains 
per layer to n + 1 chains per layer [2l[ . 

The vortex configuration for Uy < 1.04/A^ has only 
one column of vortices per pinning trough modulation, 
as illustrated in Fig. ^a). The presence of triangular 
ordering in the vortex lattice is clearly indicated by the 
structure factor S'(k), 



exp(ik ■ r,; 



(3) 



plotted in Fig. O^b). On average, the vortex- vortex in- 
teractions cancel along the x direction since the vortices 
are equally spaced a distance a apart in this direction. 
Near ny = 2.0/A^, a transition in the vortex configura- 
tion occurs in which the vortex columns begin to buckle, 
as shown in Fig. [^c). Here, the vortex- vortex interac- 
tions along the x direction no longer cancel, and at points 
in the vortex columns where a buckled vortex is adjacent 
to the column, a vortex within the column experiences 
both the force from the pinning substrate, which has 
a maximum value of Ap, as well as an additional force 
from the buckled vortex, which has a maximum value of 
Ki{a)fo. A simple estimate of the depinning force in a 
buckled region, f~f^ gives f~,, ^ Ap- Ki{a)fo = 0.395/o, 
which is close to the value of shown in Fig. [2] for 
Uy > 2.08/A^. For values of n.„ just above the buck- 
ling onset, buckling occurs in only a small number of 
places that are widely separated. Vortices near the buck- 
led locations depin at /~^, while vortices outside of the 
buckled locations do not depin until Fdc is close to Ap. 
The result is the two stage depinning process shown in 
Fig. [5] for Uy = 2.08/A^. As the vortex density is fur- 
ther increased above Uy = 2.08/A^, the distance between 
buckled regions decreases, smearing out the two step de- 
pinning process. 

For 1.91/A^ < Uy < 5.90/A^, y direction spatial cor- 
relations in the vortex configurations are lost due to the 
buckling of the vortex columns, but the pinning substrate 
maintains the spatial correlations in the x direction, re- 
sulting in smecticlike vortex structures. This is illus- 
trated in Fig. Ei;c,d) for n„ = 2.08/ X"^ and in Fig. ^e,f) 
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FIG. 5: (a) Vortex locations (black dots) and edges of the pinning troughs (black lines) in a 12A x 12A subsection of a system 
with riv — 1.04/A'^ and no applied drive. A single column of vortices is confined within each trough, (b) Corresponding 
structure factor 5(k) for n„ — 1.04/A'^. Here and in other plots of ^(k), we show a greyscale heightmap of S{kx,ky) with 
white values higher than dark values, and with the origin of coordinates at the center of the panel, is along the vertical 
axis, and ky is along the horizontal axis, (c) Vortex and pinning trough locations for a system with Uv = 2.08/A^ and no 
applied drive. In some regions, the columns of vortices have buckled and a smectic type of ordering begins to appear, (d) 
S'(k) for n„ = 2.08/A^ (e) Vortex and pinning trough locations for a system with n„ — 3.47/A^ and no applied drive. Each 
trough contains two columns of vortices, (f) S{\i) for n„ = 3.47/A'^. A double smectic ordering appears, as indicated by the 
two prominent lines in S'(k). 



for n„ = 3.47/A2. At n„ = 2.08/A2 in Fig. [^c), a 
second column of vortices begins to form in each pin- 
ning trough by means of a buckling transition, as pre- 
viously discussed. The corresponding structure factor 
in Fig. [Sl^d) shows some smectic type smearing, indicat- 
ing that the vortices are less correlated in the y-direction 
than in the x-direction. As Ny increases, the correlations 
in the y direction are further decreased, as seen in S'(k) 
for = 3.47/A^ in Fig. [5ljf). At this density, the vortex 
configuration shown in Fig. [5Je) indicates that there are 
two columns of vortices per trough, which is also reflected 
by the double band feature in S'(k). At riy = 5.21/A^, 
shown in Fig.[5]Ja), the vortex configuration still has two 
columns of vortices in each pinning trough, and the two 
band feature in S'(k) is more prominent, as illustrated 
in Fig. [HJb). As increases further, the two column 
structure is destroyed and the vortex configuration be- 
comes highly disordered in both the x and y directions, 
as seen at — 6.94/A^ in Fig. [HIc). Here S'(k) has a 
ring structure, as shown in Fig. [H^d). The effect of the 
one-dimensional substrate can still be seen in the form of 
small modulations along in Fig.[5I^d). 

At riy = 10.42/A2, plotted in Fig.EJe), a new commen- 
surate crystalline structure appears where four vortices 



fit across the pinning trough at a 45° angle so that the 
overall vortex lattice symmetry is square. This is sup- 
ported by the four prominent peaks that appear in S'(k) 
in Fig. E^f). There is a competition between minimiz- 
ing the vortex-vortex interaction energy, which favors a 
triangular lattice, and maximizing the vortex-substrate 
interaction, which favors commensurate vortex configu- 
rations. Since the difference in energy between a trian- 
gular and square vortex lattice is quite small, the com- 
mensuration effect is favored at riy — 10.42/A^, and a 
square structure appears. Since the vortex lattice is sym- 
metrical, the vortex-vortex interactions are strongly re- 
duced and the depinning transition is similar to the case 
of Uv < 1.91/A^ where /~ > /+ and the vortices depin 
more readily in the easy direction. In this case, the de- 
pinning response is elastic, but since there are multiple 
vortex columns in each trough there is a pronounced dis- 
tortion in the lattice when a dc force is applied, causing 
to be significantly lower than Ap. 

As a function of vortex density, the system thus passes 
through a sequence of triangular, smectic, disordered, 
and square vortex arrangements. These results show 
similarities to studies performed on magnetic particles 
interacting with simple symmetric periodic modulated 




FIG. 6: (a) Vortex locations (black dots) and edges of the pinning troughs (black lines) in a 12A x 12A subsection of a 
system with n„ — 5.21/A^ and no applied drive. Each trough captures two columns of vortices, (b) Corresponding S'(k) for 
rii, — 5.21/A^ showing evidence of smectic ordering, (c) Vortex and pinning trough locations for — 6.94/A^. The vortex 
configurations are disordered, (d) 5'(k) for riv = 6.94/A^ has a clear ring structure, (e) Vortex and pinning trough locations 
for Uv = 10.42/A^. The vortices form a tilted square lattice, (f) S'(k) for Uv = 10.42/A^ shows four-fold peaks. 
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FIG. 7: (Color online) Contour plot of (14) as a function 
of ac drive period t (in units of simulation time steps) and 
ac drive amplitude Fac (in units of /o) for a system with 
n„ = 0.087/A^ 



substrates, where crystalline, smectic, and partially dis- 
ordered structures were observed as a function of particle 
density [l^. The double smectic phases and the square 



lattice that we find did not appear in the magnetic parti- 
cle experiment since the relative particle densities consid- 
ered in Ref. [l9| were much lower than the relative vortex 
densities we use here. Additionally, the asymmetry in our 
potential substrate may help to stabilize certain phases 
that are not seen with simple symmetric potentials. 



IV. RATCHET AND REVERSE RATCHET 
EFFECT 

We next examine the response when an external ac 
drive is applied. When the vortex density is low, the 
vortex-vortex interactions are negligible and the system 
responds in the single particle limit. In Fig. [7] we show a 
contour plot of the average velocity (T4) obtained from 
simulations with n„ — 0.087/A^ for varied ac period r and 
ac amplitude Fac- A series of tongues where a positive 
ratchet effect occurs are clearly visible. At small r, there 
is not enough time for the vortices to respond to the ac 
drive so there is no ratchet effect. Similarly, at low Fac, 
the vortices are unable to cross the potential barrier of 
the trough, and no rectification occurs. At intermediate 
values of r and Fac, a ratchet effect can occur. In the 
first tongue, which falls at the lower left edge of Fig. [3 
the vortices move by +a in the x direction during the 
positive half of the driving period, but do not move back 
during the negative half of the period, resulting in a net 
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FIG. 8: (K) vs n„ for the system in Fig. [7] at fixed r = lO" 
simulation steps for (black circles) Fac = O.6/0 and (open 
circles) Fac = 0.4/o. Here both positive and negative ratchet 
effects occur as a function of vortex density. 



FIG. 9: (14) vs Fac for the same system as in Fig. |8] for 
varied n„. Filled squares: n„ = 1.39/A^; open squares: 
rii, = 1.74/A^; filled circles: = 2.43/A^; open circles: 
riv — 3.13/A^. The positive ratchet effect is preceded by a 
negative ratchet effect for n„ > 1.74/A^. 



rectification. On the second tongue, vortices translate 
by +2a during the first half of the period and by —a 
during the second half of the period, again giving a net 
rectification. As r and Fac increase, the ratchet effect 
on the nth tongue is characterized by motion of +na 
during the positive half period of the drive and motion 
of — (n — l)a during the negative half period of the drive. 
The tongue structure emerges since for some values of r 
and Fac, vortices undergo equal translations during both 
halves of the driving period, leading to no net motion. 
The boundary of the nth tongue where ratcheting can 
occur is given by 



nV 



Fac - F+ 



V 



1)1- 



Fa. 



F- 



(4) 



The weak pinning force F^ associated with the trough 
side of length Z+ is F^ = —Ap/2 in our case, while the 
strong pinning force F~ associated with the trough side 
of length is F~ = Ap. For the first tongue in Fig. [3 
Eq. 4 predicts ri oc {Fac + ^p/2)~^, which matches the 
simulation results. The results in Fig. [7] and Eq. 4 agree 
well with the tonguelike features found in Ref. as a 
function of t and Fac- This confirms that at low vor- 
tex densities, an ordinary ratchet effect occurs where the 
vortices move in the easy direction. 

We next consider the effect of increasing Uy. Since 
the ratchet effect in Fig. [7] is enhanced for larger peri- 
ods T, we fix the period at r = 10^ simulation steps. 
In Fig. 8 we plot (14) versus n„ for Fac — O.6/0 and 
0.4/0. For n„ < 2.60/A^, there is a positive ratchet 
effect for Fac — O.6/0 and no ratchet effect for Fac = 
0.4/0. For riy > 2.60/A^ a negative ratchet effect oc- 
curs for Fac — 0.4/0, while the positive ratchet effect for 
Fac = O-6/o decreases in size and reverses to a negative 
ratchet effect near = 3.65/A^. For both values of Fac, 



the system passes through a maximum negative value of 
(14) just below Uy = 5.21/A^. Above riy — 5.21/A^, 
{Vx) gradually increases back toward zero before cross- 
ing zero and reversing to a small positive ratchet effect 
near Uy — 10.42/A^. The features in these curves appear 
to be correlated with the vortex lattice structures pre- 
sented in Figs O and [HI At low Uy, the vortices form a 
triangular or smectic structure with only a single column 
of vortices in each pinning trough, as in Fig. [nja,b) . In 
this case, since the symmetry of the vortex configuration 
cancels or strongly reduces the vortex-vortex interaction 
force, a positive ratchet effect occurs that is determined 
only by the pinning forces of the substrate. The onset 
of the negative ratchet effect corresponds to the appear- 
ance of buckling of the columns of vortices in the pinning 
troughs, as seen in Fig. [ni^c,d). In the double smectic 
phase illustrated in Fig. [n)^e,f) and Fig. [Ua,b), the neg- 
ative ratchet effect persists. The gradual disappearance 
of the negative ratchet effect for Uy > 5.21/A^ in Fig. [H] 
is associated with the crossover from the double smectic 
state to the disordered state shown in Fig. [6l[c,d). The 
small positive ratchet effect near n„ = 10.42/A^ is corre- 
lated with the formation of the ordered square lattice as 
seen in Fig. [nie,f). 

In Fig. [9] we plot (Vx) versus Fac for a series of sim- 
ulations performed at varied Uy. For low riy < 1.74/A^, 
there is no ratchet effect until a threshold value of Fac is 
reached. (Vx) then increases to a peak value and gradu- 
ally decreases with increasing Fac- For Uy > 1.74/A^, the 
threshold value of Fac above which a ratchet effect occurs 
is lower than for n„ < 1.74/A^, and the initial ratchet ef- 
fect is in the negative direction rather than the positive 
direction. As Fac is increased, the negative ratchet effect 
reverses to a positive ratchet effect. 

Figure [TUl illustrates (Vx) in a contour plot as a func- 
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FIG. 10: (Color online) A colorscale map of (14) as a function 
of Fac and Uv For low Fac there is no ratchet effect. At higher 
Fac there is a positive ratchet effect for low reversal 
to a negative ratchet effect for intermediate riv, and another 
reversal at high n„ to a positive ratchet effect. 

tion of Fac and Uy. The positive and negative ratchet 
effect regimes are highlighted. For low Fac there is no 
ratchet effect, while for high Fac/ fa > 1.0 all ratchet 
effects are gradually reduced. The maximum positive 
ratchet effect occurs for Fac/ fo close to 1.0, the value 
of the pinning force in the strong direction. The posi- 
tive ratchet effect for Uy < 2.60/A^ does not occur un- 
less Fac/ fo > 0.5, which is the maximum pinning force 
in the weak direction. The negative ratchet effect ap- 
pears for Uy > 2.08/A^ and extends up to ny = 8.68/A^, 
with the maximum negative rectification occurring for 
3.47/A2 < < 5.21/A2 at 0.45 < Fac/fo < 0.5. For 
Uy > 9.72/A^, a reentrant positive ratchet effect ap- 
pears that is much weaker in amplitude than the positive 
ratchet effect region found for low n„. In general, the 
ratchet effect in either direction is reduced for increas- 
ing Uy since the effects of the substrate are gradually 
washed out by the increasing vortex-vortex interaction 
forces. This result shows that vortices interacting with 
simple asymmetric pinning potentials can exhibit multi- 
ple ratchet reversals similar to those seen for the more 
complex 2D pinning array geometries [l3|. It is likely 
that for higher Uy , higher order reversals will occur when 
the vortices again form a smectic phase with five or more 
vortex chains per pinning trough. We are not able to ac- 
cess the high vortex densities where these reversals would 
be expected to appear. The higher order ratchet effects 
would likely continue to decrease in amplitude with in- 
creasing field. 

The positive or regular ratchet effect that occurs when 
the vortices move in the easy direction is straightforward 
to understand at low fields based on a single particle pic- 



(b) (d) 

FIG. 11: (a) Vortex locations (black dots) and edges of the 
pinning troughs (black lines) in a 12A x 12A subsection of a 
system with — 2.08/A^ at Fac = 0.5/o when the nega- 
tive ratchet effect occurs. The positive portion of the driving 
period, with F^^' — UFac, is shown. The vortices form two 
columns in every trough, (b) The corresponding structure fac- 
tor ^(k) for positive drive. The overall disorder in the vortex 
configuration is indicated by the ring structure, (c) Vortex 
and pinning trough locations for the same system during the 
negative portion of the driving period, with Ff^* — —UFac- 
The vortices are arranged in buckled columns, (d) The cor- 
responding ^(k) for negative drive shows a smectic type of 
ordering. 



ture. The first reversal of the ratchet effect for increasing 
Ny is related to the vortex-vortex interactions. This is 
more clearly seen from Fig. 1111 where we plot the vor- 
tex positions and corresponding structure factor during 
both halves of the driving cycle for n„ = 2.08/A^ and 
Fac = 0.5/o where a negative ratchet effect occurs. Dur- 
ing the positive drive portion of the ac cycle, the vortices 
are roughly evenly spaced in a disordered arrangement, 
as shown in Fig. [TlT a.b). When the negative drive is 
applied, the vortices form a clear smectic structure with 
roughly one vortex column per pinning trough, as illus- 
trated in Fig. [TlT c.d). In this case, the vortex configura- 
tion is asymmetric due to the appearance of buckling in 
the vortex columns. At a buckled location, an extra vor- 
tex appears on the positive x side of the vortex column 
(to the right of the column) in the same trough. These 
extra vortices effectively push the adjacent vortices in the 
column over the potential barrier in the hard or negative 
X direction. 

If the extra vortex is located a distance a from the 
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column, then a portion of the vortices in the column can 
escape over the potential barrier in the hard direction 
when the following condition is met: 

Fac + Ki{a)fo> Ap. (5) 

When Ap/ fo — 1.0 and Ki{a) = 0.6, a negative ratchet 
effect is predicted to occur when Fac/fo > 0.4, which 
is what is observed in Fig. [HI It could also be argued 
that when the drive is in the positive direction, the ex- 
tra vortices should be pushed by the column of vortices 
and should more easily move in the positive direction, 
allowing the conditions of Eq. 5 to be met in the posi- 
tive direction as well. In Fig. [TlT a) we show that this 
does not occur because the vortex arrangement is quite 
different during the positive drive portion of the cycle. 
The vortex-vortex interaction term is not present since 
the vortices do not form the smectic structure that is 
seen when the driving is in the negative direction. In- 
stead, the vortices form two columns within each pinning 
trough, allowing the individual vortices to remain evenly 
separated in both the x and y directions, and prevent- 
ing an individual vortex from pushing another vortex out 
of the trough. The vortex structure for positive drive is 
more symmetrical, as seen in S'(k) in Fig. IllT b) which 
has a smeared crystalline or disordered structure. This 
symmetry reduces the vortex-vortex interaction forces. 

A ratchet effect in the positive direction occurs when 
the following equation is satisfied: 

Fac > Ap/2. (6) 

If Fac/fo = 0.45 and Ap/fo = 1.0, then the conditions 
of Eq. 5 are met but the conditions of Eq. 6 are not, so 
the negative ratchet effect can occur but not the positive 
ratchet effect. Once Fad fo > 0.5, it is also possible for 
the positive ratchet effect to occur and the two ratchet 
mechanisms will compete, decreasing the magnitude of 
the negative ratchet effect as seen in the phase diagram 
of Fig. [in] for higher Fac- As n„ is increased, the smectic 
structure is gradually lost as shown in Fig. [S]^c,d), so the 
vortex-vortex interactions required to produce the nega- 
tive ratchet effect are reduced and the negative ratchet 
effect is diminished. Once the vortex lattice is again sym- 
metrical, as in Fig. [ni^e,f), the positive ratchet effect ap- 
pears. This result suggests that in regions where a re- 
versed ratchet effect appears, the vortex lattice structure 
must exhibit some asymmetry when driven in the hard 
direction for the reverse ratchet effect to occur. 

We have also investigated the reversal of the ratchet 
effect as a function of the pinning strength Ap for fixed 
= 3.13/A2 and r = 10"*. In Fig. [H we plot (t4}iV„ 
versus Ap for Fac ranging from 0.4/o to O.S/q. For all 
values of Fac^ as Ap increases a positive ratchet effect 
appears. This is because the stronger Ap forces the vor- 
tices to form completely ID unbuckled columns inside 
the pinning troughs. For lower Ap, the vortex columns 
can buckle, providing the unbalanced vortex- vortex inter- 
actions required to produce the negative ratchet effect. 
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FIG. 12: {Vx)N^ vs Ap for the same system as in Fig. [TO] 
with nu = 3.13/A^. Squares: Fac = 0.3/o; circles: Fac = 
0.4/o; up triangles: Fac = 0.5/o; down triangles: Fac = O.6/0; 
diamonds: Fac ~ 0.7/o; stars: Fac = O.8/0. For large Ap 
only a positive vortex ratchet effect occurs as the vortices 
form a single column in each pinning trough. As the pinning 
strength is reduced, buckling of the columns begins to occur 
and a negative ratchet effect appears. 

The negative rectification is enhanced for both low Ap 
and low Fac- Changing the strength of Ap in a single 
superconducting sample is difficult; however, a geometry 
similar to that shown in Fig. [1] could be created using 
tunable optical traps for colloidal particles driven with 
an ac electric field. 



V. SUMMARY 

In conclusion, we have shown that the vortex ratchet 
system proposed by Lee et al. in Ref. 5] for vortices in- 
teracting with simple periodic asymmetric pinning sub- 
strates can exhibit a series of ratchet reversals similar to 
those observed for 2D asymmetric periodic pinning ar- 
rays. Unlike most vortex ratchet geometries considered 
to date, the ratchet effects presented here can only be ex- 
plained by a 2D description of the vortex configuration, 
and not by a simple ID model. As a function of vortex 
density, a rich variety of vortex phases can be realized in 
this system including triangular, smectic, disordered, and 
square lattices. The dc critical depinning forces exhibit a 
reversible diode effect as a function of vortex density. At 
low magnetic fields, the depinning force is higher in the 
hard direction, corresponding to the regular diode effect, 
while at higher magnetic fields, the depinning force is 
lower in the hard direction. Another reversal of the diode 
effect occurs at even higher fields. The reversed diode 
effect appears when the symmetry of the vortex config- 
uration changes for different directions of dc drive, as in 
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the case of the smectic and disordered phases. Here, the 
importance of vortex-vortex interactions differs for the 
two drive directions. At densities where the vortex lat- 
tice is highly ordered, a regular diode effect occurs where 
the vortex depinning force is lower in the easy direction. 
The vortex lattice structures are also directly correlated 
with the sign of the ratchet effect that occurs under an 
ac drive. The highly symmetric structures ratchet in the 
easy direction since the vortex-vortex interactions effec- 
tively cancel, causing the sign of the ratchet effect to be 
determined only by the asymmetry of the pinning sub- 
strate. In the smectic states, the vortex-vortex interac- 
tions become relevant and can produce a negative ratchet 
effect, which can be predicted with force balance argu- 
ments. We have shown the dependence of the reversible 



ratchet effect on ac amplitude, vortex density, and pin- 
ning strength. This type of ratchet could be also real- 
ized for other kinds of systems composed of collections of 
repulsively interacting particles and simple asymmetric 
substrates, such as colloid systems. 
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